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In this paper we study an anisotropic model generated from a Bianchi III metric, generalization of 
Godel's metric, which is an exact solution of Einstein's field equations. In particular, we analyse 
supernova la data, namely the SDSS sample calibrated with the MLCS2k2 fitter, verifying in which 
limits of distances and redshifts the anisotropy of the model could be observed, and in which limits 
the model approaches the ACDM model. We found that redshifts above z « 2 has a particular 
importance in setting the point at which the anisotropy would be noticed, as well as the point at 
which the present model begins to diverge from the ACDM. We conclude that data above this limit, 
as well as an increasing accuracy in measuring these distances, might indicate the existence of such 
an anisotropy in the universe. 

PACS numbers: 



The search for evidence of possible anisotropies in the universe has been studied by several authors [TJ HJ |31 HI [5] . 
The importance of these studies lies in the fact that the discovery of a possible anisotropy could lead to a better 
understanding of the mechanism of the Big Bang and of the small temperature fluctuations detected in the cosmic 
microwave background radiation (CMB) !6 : . In addition, such a discovery would bring down the belief in the validity 
of the cosmological principle that is the basis of many interpretations of observations. 

In order to detect possible signatures of anisotropy, in this paper we consider an anisotropic model constructed 
from an anisotropic (but homogeneous) metric with expansion, which is an exact solution of Einstein's equations. 
This metric, whose first particular case was proposed by Godel in 1949 [7], is classified as Bianchi III and was studied 
in a more general form by M. Rebougas and J. Tiomno jS], and by V. Korotkii and Y. Obukhov [5] (from now on we 
will call it the RTKO metric) . In general it has non-zero rotation (but zero shear) , as well as a conformal expansion. 
Here we will take the particular case where the rotation is made zero [TOl HI] • I n this case, it becomes 1 



where r\ is the conformal time. 

As is well known, a perfect fluid does not lead to anisotropies. Therefore, in order to account for them, a scalar field 
is introduced, which will be responsible for creating a pressure difference in the preferred direction, which we define as 
being the z direction. In addition, we will introduce a cosmological constant responsible for the present acceleration 
in the expansion of the universe. Although the metric is anisotropic, it provides an strictly isotropic CMB (at the 
background level) due to the existence of a conformal Killing vector parallel to the fluid 4- velocity [9l [12] . In what 
concerns the distance-redshift relations, we will consider here 288 supernovae la with redshifts between 0.0218 and 
1.551, collected by the Sloan Digital Sky Survey II (SSDS-II) [T3] and calibrated with the Multicolor Light Curve 
Shapes (MLCS2k2) fitter [H]. 

The work is organized as follows. In section 2 we present the main properties of the RTKO metric. In Section 3 we 
build the correspondent cosmological model. In Section 4 we show the method of determining the model parameters. 
In section 5 we obtain the results for the anisotropic model and, for comparison, for the ACDM model. Finally, in 
section 6 we present our conclusions. 
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II. THE RTKO METRIC 

Before proceeding to the construction of the anisotropic model, let us discuss some properties of the RTKO metric. 
If we apply the coordinate transformations 

e x — cosh r + cos <p sinh r, (2) 
ye x = sin ip sinh r, (3) 

to metric ([I]), we put it in the cylindrical form 

ds 2 = a 2 (r])(dr) 2 - dr 2 - sinh 2 rdip 2 - dz 2 ), (4) 

which is invariant under rotations around the z-axis. Now, making use of the transformations 

r = x sm #, (5) 
z = x cos 8, (6) 

the above line element can be expressed in spherical coordinates as 

ds 2 = a 2 (7])[dri 2 - d X 2 - X 2 d6 2 - sinh 2 ( X sin 6) dip 2 ). (7) 

The anisotropy is clear in the last term. It will lead to an angular dependence in the angular-diameter distance, 
as we shall see. Interestingly enough, the metric expressed in the later form reduces to the spatially flat Friedmann- 
Lemaitre-Robertson- Walker (FLRW) space-time in the limit of small distances. Indeed, taking the expansion of the 
hyperbolic sine function up to first-order terms, we obtain 

ds 2 a 2 (r])[dri 2 - d X 2 - X 2 {d9 2 - sinh 2 6dip 2 )}, (8) 

i.e., the flat FLRW metric expressed in spherical coordinates. As a consequence, the anisotropic model to be built 
will reduce to the standard model in the limit of small distances. 
The RTKO metric has Killing vectors 

=9 X - yd v , £ (2 ) = d y , £ (3) = d z , (9) 
which classify it as a spatially homogeneous Bianchi type III metric. It also has a conformal Killing vector 

e conf - # . (io) 

The 4-velocity of a comoving fluid in this space-time is given by 



dr^ ^ £ M t 
ds a(rj) a(r]) 

Let us consider now the Einstein field equations 



(11) 



K - \KR - 1?. (12) 

When we use metric ([!]), we obtain for the energy- momentum tensor the diagonal components 

T V = 3a' 2 -a 2 , (13) 
Tla 4 = Tia 4 = 2aa" -a' 2 , (14) 

r 3 



Tfa 2 = T\a 2 - 1, (15) 



while all the non-diagonal ones are zero. Here, the prime denotes derivative with respect to the conformal time. In 
the next section we discuss the cosmological models arising from the above metric. 
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III. THE ANISOTROPIC MODEL 



Let us define 



T-lU 

-t n 



e T? 



-pi 



(16) 



as the total energy density and pressures of the universe content, where the index i takes the values 1, 2 and 3, 
representing the pressures in the directions x,y e z, respectively. Equations ( 13 )-( L5 1 then become 



4 o '2 2 

ea = 6a — a , 



4 4 /2 r\ // 



P3a 



Pi a 



1. 



(17) 
(18) 
(19) 



This set of equations shows clearly the anisotropy of the metric, since the pressure along direction z differs from 
the others two by an additional term 1/a 2 . The usual components of the universe (radiation, matter and cosmological 
constant) can not generate such an anisotropy. Let us show that it can be generated, in a self-consistent way, by a 
massless scalar held minimally coupled to gravity. Such a scalar held satisfies the Klein-Gordon equation 



1 



(20) 



and has an energy-momentum tensor given by 



(21) 



Here g is the determinant of the metric tensor g^, and 8tf is the delta of Kronecker. The commas and semicolons 
refer, respectively, to ordinary and covariant derivatives. 
Our scalar held will be given by 



$(z) = Cz, 



(22) 



where C is a constant. It is easy to verify that it is solution of equation (20 1. On the other hand, the non-zero 
components of its energy-momentum tensor (the diagonal ones) are 



T ° = —Ty = -Tf = T| = 



C*_ 
2a 2 ' 



which lead, through (16), to 



P W = 



C*_ 
2a 2 ' 



(23) 



(24) 



(«) (») 
Pi = V\ 



(«) 

-P3 



2a 2 ' 



(25) 



where the upper index (s) refers to the contribution of the scalar held. We note that the pressures are anisotropic, 
and the energy density falls with the square of the scale factor. Using this fact we can work out the anisotropy present 
in Einstein's equations. 
Let us dehne 



e = e + e (s) , Pi = Pi + pf\ 



(26) 



where the bar refers to the isotropic content (radiation, matter and cosmological constant). Equation (15) takes the 
form 



Pz a 2 = Pl a 2 — C 2 + 1. 



(27) 



We then see that the anisotropy of the model is absorbed by the scalar held if (and only if) C 2 = 1. In this case we 
have pj = p" 2 = p~ 3 = p, and the scalar held is simply given by 



$0) = ±z. 



(28) 
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The energy density and pressures are now given by 



2a 

and the Einstein equations can be rewritten as 



= d?> (29) 

pW= p W = -pW = _-L ! (30) 



ea 4 = 3a' 2 -^a 2 , (31) 
pa 4 = a' 2 -2aa"+^. (32) 



Notably, these equations are precisely those of an open FLRW model with curvature k = —1/2. In fact, from (31 ) 
we obtain the Friedmann equation 

where we have defined the Hubble parameter H = a' /a 2 . 
On the other hand, from the continuity equation 

- + 'SH(e+p) = (34) 
a 

we can write the energy densities of radiation, non-relativistic matter and the cosmological constant respectively as 
e r = -4r, e rn — and e\ — A, with A, B and A constants. The total energy density of the isotropic content is then 

A " A. (35) 

a 1 * a J 



The Friedmann equation ( 33 ) assumes now the form 

, A B A 1 
^=3^ + 3V + 3 + 2^ ^ 

It can also be written as 

f} (s) = l-tl r -Q m - A , (37) 
where, for matter, radiation and the cosmological constant, we have defined flj = ej/'SH 2 (j = r, m or A), and 

QM s 2ik?- < 38 > 

Owing to its particular dependence on the scale factor, we will call this last term the "curvature" of the anisotropic 
model. However, let us point out that it contains, besides the space-time curvature itself, the contribution of the 
energy density of the scalar field responsible for the anisotropy. 
By using the relation 

1 + z-^, (39) 



we obtain, from ( 36 1 and ( 37 ) , 
H(z) 



= E(Z) = vtoro(l + *) 4 + + Zf + (1 - tiro - fi m0 - fi A )(l + zf + fi A) (40) 

where the index refers to the present time. On the other hand, the deceleration parameter is defined by 



or, as a function of the redshift, by 



while the age parameter is given by 



a a 

q=\--p>, (41) 



fl + z) dH , 

^)=- i+ W^' (42) 



H to= / — 777 -,dz. (43) 
Jo l + zH(z) 
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IV. PARAMETERS ESTIMATION 



In order to establish a comparison between the anisotropic model and the standard one (here denominated fc-ACDM 
to include the spatially curved cases), we have performed an analysis of 288 supernovas la compiled by the SDSS 
Supernova Survey with the fitter MLCS2k2. The best-fit values of the free parameters of both models were determined 
by means of the usual x 2 statistics. 



A. The fc-ACDM case 



In this case we use the modulus distance 

5 log 

/jl(z) — <* 5 log 



AO 



v4^ 



AO 



dz 



51og(» +42.38 (fc = +l) 



(1 + *) 



dz 

W) 



5 log 



\ _i_ z / 

/its 1 sinh VPfcOi / - r 



o E(z), j 

- 51og(/i) + 42.38 (fc = 0) 

dz 



- 51og(» +42.38 (fc = -1), 



where h = ffo/100 km/s-Mpc, and 

e{z) = v/n m0 (i + ^) 3 + (i - ft r0 - n m0 - «a)(i + z? + Pi. 

After marginalising /i and fixing Sl r o with the CMB constraints, the free parameters are £l m o and Pa- 



(44) 



(45) 



B. The anisotropic case 

The FLRW metric is isotropic. This means that, given two galaxies with the same redshift, they are at the same 
luminosity distance from us. In the case of an anisotropic metric that is not true. For a given redshift we can obtain 
different distances, i.e., the brightness depends on the distance and on the angle of observation. We need therefore to 
obtain the correct luminosity distance for our metric. 

With this goal we will define the angular-diameter distance by [4[ [T2] 



2 dA p 

dA =ln- 



(46) 



where dA p is the proper-area element formed by the solid angle dft = smOdOdip. The reason for this choice is the fact 
that, in an anisotropic metric, the usual definition 



d A 



dl_ 
da 



(47) 



(where da is the angle associated to the distance dl between two observed points) would lead to different distances dl 
for the same d A as the observation angle 9 changes. 

The area element in the anisotropic case can be obtained by taking our spatial metric and fixing Xi which leads to 
the line element 



dl 2 = a 2 (T])[x 2 de 2 + sinh 2 (x sin 6») V] 
from which we can see that the proper-area element is given by 

dA p {rj) — a 2 (ry)xsinh(x sin 9) d6 'd<p, 

which, after dividing by <if2, leads to 

sinh(xsin6')" 1 ly ~ 



d A {rj) = a(v)x 



Xsin9 



(48) 
(49) 

(50) 
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In order to perform an observational analysis, we can rewrite equation (501 in terms of the redshift, 

sinri(xsinfl)" 1 11 ' 



d A (z) = 



aoX 
l + z 



Xsmt 



6e (0,tt). 



On the other hand, the coordinate \ can be obtained by taking a null geodesic, 



dx = drj =4> x = 



1 



where we have defined 



Z[z) 



a Q H ( 



dz' 



■Z(z), 



o E(z') 



Therefore, the luminosity distance will be given by 

(l+z)Z(z) 



d L (z) = (i + z yd A (z) = 



H 



v /|r2 fe0 |Z(z)sin( 



-,1/2 



G (0,tt). 



(51) 



(52) 



(53) 



(54) 



We can note that, in the limit 9 — > or x "C 1, the luminosity distance reduces to that of the spatially flat FLRW 
model. Let us also remind that |f2fco| = (see (38)), which permits to rewrite the above equation in terms of Qq S ^ . 

In order to fix the free parameters of the present model, |f2 m o, fi^ | (after marginalizing h), we use the angular 

average of the luminosity distance in the interval [0,7r], for each value of z. The average modulus distance will be 
given by 



■Jl(z) = 5 \og[H d L (z)} - 5 \og{h) + 42.38, 



(55) 



where di,(z) is the average of (54) in the interval 9 € [0,7r] 



V. RESULTS AND DISCUSSION 



With the procedure described above, we have obtained the best fit for the fc-ACDM model with parameters £l m0 — 
0.39, 11a = 0.58 and flko = 0.03, with Xr — 0.84. Therefore, in the absence of other priors, the data favors a slightly 
negative curvature. Within 2er confidence level, we obtain 0.16 < 51 m o < 0.58 and 0.12 < fl\ < 0.95. For the 
curvature we have, within the same confidence level, —0.51 < f^o < 0.71. 

For the anisotropic case, the best fit corresponds to S7 Tn0 = 0.38, £1a = 0.58 and fi^ = 0.04, with Xr = 0.84. The 
2(7 intervals are < n m0 < 0.48, 0.04 < n A < 0.67 and < < 0.95. The confidence regions for both models 
are presented in Fig. 1. In Table I we compare the parameters best values for both models, as well as the present 
deceleration parameters qo and the age parameters Hoto, obtained from equations (42) and (43), respectively. 

In order to set the level of agreement between the anisotropic model and the observational data, in Fig. 2a we show 
a comparison between the observed distance moduli and the theoretical predictions. The solid lines correspond to 
the maximum {9 = 7r/2), average and minimum (9 = 0) values of the distance modulus. In this comparison we have 
used Q( s ' — 0.52, which is the maximum "curvature" allowed (within the 2a level) when we fix the matter density 
parameter in 51„ l0 = 0.25. We can see that, for the range of redshifts available in SDSS, the theoretical predictions 
are within the observed limits. 

As shown above, the best fit requires a small curvature (a spatially open universe) in the fc-ACDM model. There 
is then space for confusing a possible anisotropy with an actual space curvature. In order to verify in which range of 
redshifts the two models may be confused, we show in Fig. 2b the average distance modulus as a function of z in both 



Model 






«a 


x, 2 


Qo 


Hoto 


fc-ACDM 


0.39 


0.03 


0.58 


0.840 


-0.386 


0.888 


Anisotropic 


0.38 


0.04 


0.58 


0.840 


-0.384 


0.890 



TABLE I: Parameters best values for the fc-ACDM and anisotropic models fitted with the sample SDSS (MLCS2k2). 
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FIG. 1: Confidence levels la, 2a and 3a for the fc-ACDM (white) and anisotropic (blue and white) models with SDSS 

(MLCS2k2). 
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(a) (b) 

FIG. 2: (a) Distance modulus versus redshift for the anisotropic model (minimum, medium and maximum) and (b) for the 

models with best fit parameters for z < 2. 

models, for the corresponding best-fit values of their parameters. We can note that both agree with the observational 
data, and that the curves superpose even for redshifts above z ~ 2 (see Fig. 3a). However, when we use the maximum 
curvature allowed (within 2a level) for f2 m o = 0.25 (which corresponds to — 0.52 and flko = 0.64), the models 
predictions start to diverge (see Fig. 3b), indicating that, above the upper limit of the used data, the test is more 
sensitive to a possible anisotropy, which could, in principle, be detected. 

VI. CONCLUSION 

We have studied the observational viability of an anisotropic cosmology based on a Bianchi III metric which 
generalizes the Godel solution, in the particular case where the cosmic rotation is made zero. As shown elsewhere, 
such a metric leads to an exactly isotropic CMB (at the background level) and, in addition, it is an exact solution of 
the Einstein equations in the presence of an anisotropic, minimally coupled scalar field. 

This metric induces an anisotropy in the luminosity distance of galaxies with same redshifts but observed in different 
directions. This effect can, in principle, be used to detect a preferential axis in the sky through the observation of 
supernovas la at high redshifts. In the present work we have made use of the SDSS supernova sample calibrated with 
the MLCS2k2 fitter, since it is less dependent on a ACDM fiducial model as compared to SALT or SALT II fitters. 
For comparison, we have also analysed the ACDM model, with any possible curvature. 

Our analysis have shown that, in both cases, a small negative curvature (or anisotropy) is allowed by observation. 




On the other hand, the present data do not distinguish between the two models, even when we consider the maximum 
permitted curvature/anisotropy. A more robust conclusion about the existence of an anisotropy at cosmological scales 
may, however, be reached when supernovas data with higher redshifts are available. 
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